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Abstract. We study the limiting behavior of a singularly perturbed Schrodinger- 
Poisson system describing a 3-dimensional electron gas strongly confined in the 
vicinity of a plane (x, y) and subject to a strong uniform magnetic field in the 
0^ ' plane of the gas. The coupled effects of the confinement and of the magnetic 

c 2 ^ , field induce fast oscillations in time that need to be averaged out. We obtain 

' at the limit a system of 2-dimensional Schrodinger equations in the plane {x,y) 

(N ■ coupled through an effective selfconsistent electrical potential. In the direction 

5— ( , perpendicular to the magnetic field, the electron mass is modified by the field, as 

Oh' the result of an averaging of the cyclotron motion. The main tools of the analysis 

I are the adaptation of the second order long-time averaging theory of ODEs to 

. our PDEs context, and the use of a Sobolev scale adapted to the confinement 

' operator. 

<. 

r-| '. 1. Introduction 

» . 

j3 ' l-l- The singularly perturbed problem. Many electronic devices are based on 

the quantum transport of a bidimensional electron gas (2DEG) artificially confined 
in heterostructures at nanometer scales, see e.g. [21 HI EOl EI] • In this article, we 
CN ■ derive an asymptotic model for the quantum transport of a 2DEG subject to a 

^ i strong uniform magnetic field which is parallel to the plane of the gas. The aim of 

this paper is to understand how the cyclotron motion competes with the effects of 
O . the potential confining the electrons and the nonlinear effects of the selfconsistent 

Poisson potential. Our tool is an asymptotic analysis from a singularly perturbed 
CN ' Schrodinger-Poisson system towards a reduced model of bidimensional quantum 

, transport. In particular, we generalize in this context the notion of cyclotron ef- 

fective mass, usually explicitely calculated in the simplified situation of a harmonic 
confinement potential [20l [28] . 

Our starting model is thus the 3D Schrodinger-Poisson system, singularly per- 
^ ' turbed by a confinement potential and the strong magnetic field. The three- 

dimensional space variables are denoted by (x, y, z) and the associated canonical 
basis of is denoted by {ex, ey, Cz)- The particles are subject to three effects: a 
confinement potential depending on the z variable, a uniform magnetic field applied 
to the gas along the Cy axis, and the selfconsistent Poisson potential. Given a small 
parameter e > 0, which is the typical extension of the 2DEG in the z direction, our 
starting model is the following dimensionless Schrodinger-Poisson system: 

idt"^^ = ^ [-dl + B'^z^ + Vc{z)) - ^2iBzdx^^ - /^x,y^^ + V'^' , (1.1) 

^%^,x,y,z) =^Q{x,y,z), (1.2) 

V%t,x,y,z) = ^^.m\ (1.3) 

where we have denoted 



in 
o 



0^ 
O 



r'^{x, y,z) = ^/x"^ + y"^ + e^z"^. (1.4) 
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The scaling is discussed in the next subsection. This system describes the transport 
of electrons under the action of: 

- The applied confinement potential -^Vciz), nonnegative, such that Vc{z) — > 
+00 as \z\ +00. The precise assumptions of this potential are made below 
in Assumptions 11.11 and 11.21 



- The applied uniform magnetic field ^Cy (with B > fixed), which derives 
from the magnetic potential ^BzCx- We have chosen to work in the Landau 
gauge. 

- The Poisson selfconsistent potential . 
Note that (jl.ip is equivalent to 

idt^^ = ^{-d'^^ + v^iz))^' + (^id^-^^ -d^^' + v'^' . (1.5) 

The goal of this work is to exhibit an asymptotic system for (jl.ip . (|1.2p . p.3p as 
e ^ 0. 

Let us end this subsection with short bibliographical notes. In a linear setting, 
quantum motion constraint on a manifold has been studied for a long time by 
several authors, see [ISl HU EH [30] and references therein. Nonlinear situations were 
studied more recently. The approximation of the Schrodinger-Poisson system with 
no magnetic field was studied when the electron gas is constraint in the vicinity 
of a plane in [3, [25] and when the gas is constraint on a line in [5]. When the 
nonlinearity depends locally on the density, as for the Gross-Pitaevskii equation, 
asymptotic models for confined quantum systems were studied in [H [Gj [T^. In 
classical setting, collisional models in situations of strong confinement have been 
studied in |T7] . Finally, let us draw a parallel with the problem of homogenization 
of the Schrodinger equation in a large periodic potential, studied in [T] and f^H]. 
At the limit e — > 0, as noted above, we will obtain an homogenized system which 
takes the form of bidimensional Schrodinger equations with an effective mass in 
the X direction. However, this phenomenon is due to an averaging of the cyclotron 
motion induced by a strong magnetic field, and is not exactly the same notion as 
the usual effective mass for the transport in a lattice or in a crystal. Nevertheless, 
it is interesting to observe that the scaling used in [TJ [29] in the case of a strong 
periodic potential is similar to the strong confinement scaling used in the present 
paper. 

1.2. The physical scaling. In order to clarify the physical assumptions under- 
lying our singularly perturbed system, let us derive (11.11) . (11.21) . (|1.3I) from the 
Schrodinger-Poisson system written in physical variables. This system reads as 
follows: 

1 / . „ eB ^ ^ 



ihdt^ = [ihV zcx * + eVc*-FeV*, (1.6) 

2111 \ c J 

V = - . / ^ (1.7) 

Each dimensionless quantity in p.ip . p.2p . p.3p is the associated physical quantity 
normalized by a typical scale: 

X = ^, y = ^, . = ^, 1^1^ = ^, V. = l^,V^ = l,B = 2. (1.8) 

x' ^ y' z' ' ' N Vc V B ^ ' 

Now we introduce two energy scales in this problem: a strong energy Econf, which 
will be the energy of the confinement in z and of the magnetic effects, and a transport 
energy Etransp, which will be the typical energy of the longitudinal transport in (x, y) 
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and also of the selfconsistent effects. We introduce the following small dimensionless 
parameter: 

Econf ) 

Then our scaling assumptions are the following. We set to the scale Econf the 
confinement potential, the magnetic energy and the kinetic energy along z: 

Econf .= eVc = -m { — ] z^ = :—2 1-10 
2 \mc J 2mz'^ 

and we set to the scale Etransp the selfconsistent potential energy, the kinetic energies 
along X and y and we finally choose a time scale adapted to this energy: 

e 2mx^ 2my^ t 

By inserting ([LEI) in dTH)), ([r7|, then by using ([TH]), (fTTOD and (fTTTl) . we obtain 
directly our singularly pertubed problem p.lj) . p.3p . Note that p.lOjl and (II. lip 
imply that e is also the ratio between the transversal and the longitudinal space 

scales: _ _ 

z z 

X y' 

1.3. Heuristics in a simplified case. In this section, we analyze a very simplified 
situation where analytic calculations can be directly done. We assume here that Vc 
is a harmonic confinement potential and we neglect the Poisson potential V^. We 
formally analyze the heuristics in this simplified case, that will be further compared 
to our result obtained in the general case. 

We thus consider here a new system, similar to (II. ip where we prescribe Vc{z) = 
oP'z^ ^ a > and where the Poisson potential is replaced by 0: 

idt^^ = ^ {-dl + (a^ + b'^)z'^) - -2iBzd^^' - A^,j,^= , (1.12) 

^"(0,x,y,z) =^o(x,y,^). (1.13) 
In this situation, there is a trick which enables to transform the equation. Indeed, 
by remarking that 

-dl + (a^ + B'^)z^ - 2iBezd^ - e^dl 

2 • 

Q,2 + B2 —J „2 + ' 

we obtain that p.l2|) is equivalent to 

52 , / ,2 , d2n " 



-dl + {a^ + B^)[z--,^^^ied.A --^e^dl 



-di + {a' + B^)[z-^^^ted, 



a' 



2 _|_ ^2 a; - 



(1.14) 

Introduce now the following operator: for a function u E L^(M^), we set 

(e^n)(x,y,2;) = (^^xu{C,y,z + -^^^^-^eC)^ , 

where Tx denotes the Fourier transform in the x variable. Note that this operator 
0^ is unitary on L^(M^) and commutes with dx and dy. Hence, we deduce from 
(|1.14|) and by direct calculations that the function u'^ = B^^''^ satisfies the following 
system: 

idtu^ = \h,u' - - d^u', u'{t = 0) = O'^q , 
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where _ 

Let us now filter out the oscillations by introducing the new unknown 

= exp(it5^/e2)n^ 

Again, the operator exp{itHz/e'^) commutes with dx, dy and, finally, the following 
equation is equivalent to (|1.12l) : 

2 

= - 9y^"^ ^'(* = 0) = ©'^0 • (1.15) 

Q.Z _|_ y 

As e ^ 0, it is not difhcult to see that, for sufficiently smooth initial data, we have 
B'^^o — > ^0- Therefore, one can show that, in adapted functional spaces, we have 
— > $ as e ^ 0, with ^ solution of the limit system: 

idt^ = -^^^^dl^-d^y^, ^t = 0)=^0. (1.16) 

This equation is a bidimensional Schrodinger equation with an anisotropic operator 
that can be interpreted as follows. Whereas, as expected, the dynamics in the y is 
not perturbed by the magnetic field (since it is parallel to y), in the x direction the 
electrons are transported as if their mass was augmented by a factor " > 1. 
This coefficient is called the (dimensionless) electron cyclotron mass [20l 

In this article, the model that we want to treat is the nonlinear system (II. 
p.3p . with a general confinement potential Vc instead of a^z^ and the selfconsis- 
tent Poisson potential. Consequently, it is not possible to simplify the equation 
p.ip by the above trick. Moreover, the potential depends on the z variable 
and on the function itself. Therefore, one has to be careful for instance when 
filtering out the fast oscillations by applying the operator exp{itHz/e'^), since in 
this nonlinear framework some interference effects between the elementary waves 
might appear. In this article, we present a general strategy that enables to over- 
come these difficulties. The strategy will be inspired from [6] where the nonlinear 
Schrodinger equation under strong partial confinement was analyzed. Two main 
differences appear here. First, the Poisson nonlinearity is nonlocal, which requires 
specific estimates. Observe that, at the limit e — > 0, the nonlinearity in the present 
paper reads * JlV'pdz and does not depend on z. This makes an important 
difference with the case of |6j, in particular no resonance effects due to the nonlin- 
earity will appear. Second, the magnetic field induces in (II. ip a singular term at 
an intermediate scale ^ between the confinement operator (at the scale -ij) and the 
nonlinearity (at the scale ^). Hence, compared to [6], the average techniques have 
to be pushed to the order two and resonance effects will finally appear here due to 
this magnetic term. 

1.4. Main result. Consider the system (11.11) . (11.2(1 . (II. 3p . We assume that the con- 
finement potential Vc satisfies two assumptions. The first one concerns the behavior 
of this function at the infinity. 

Assumption 1.1. The potential Vc is a C°° nonnegative even function such that 

c?\z^ < Vc{z) < C\z\^'^ for \z\ > 1, (1.17) 
where a > 0, M > 0, and 

\BM^)\^o(\.r"'), MlM = o(i)/„,„,ueN-, (1.18) 

Vc{z) V / Vc{z) 

as \z\ +00, where M' > 0. 
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Note that a smooth even potential of the form Vc{z) = C\z\^ for \z\ > |zo|, with 
C > 0, s > 2, satisfies these assumptions. In particular the harmonic potential 
Vc = a^z^ fits these conditions. 

Let us discuss on the assumptions. The assumption that the function Vc{z) is even 
is important in our analysis, see e.g. Step 4 in subsection 11.51 The left inequality in 
the first condition (jl.lTp implies that Vc tends to +oo as |z| — > +oo. The fact that 
VciA > c?z^ is not essential in our analysis but simplifies it (see below, it allows to 
give a simple characterization of the energy space related to our system). As it is 
well-known [26], the spectrum of operator B.^ defined by 

H, = -dl^B^z^ ^Vc{z). (1.19) 

is discrete, when Kz is considered as a linear, unbounded operator over L^(R), with 
domain 

D{Hz) = {u G l2(M), G l2(M)}. 

The complete sequence of eigenvalues of Hz will be denoted by (£'p)pgN, taken 
strictly increasing with p (recall indeed that in dimension 1 the eigenvalues are sim- 
ple), and the associated Hilbert basis of real- valued eigenfunctions will be denoted 
by (xp(-z))peN- The right inequality in (jl.lTp and the second condition p.l8|) are 
more technical and are here to simplify the use of a Sobolev scale based on the op- 
erator Hz^ which is well adapted to our problem. More precisely, these assumptions 
are used in Lemma [231 

The second assumption on concerns the spectrum of the confinement operator 
Hz. 

Assumption 1.2. The eigenvalues of the operator- Hz defined by (I1.19P satisfy the 
following property: there exists C > and riQ G N such that 

VpGN, Sp+i-Sp>C7(l+p)-"o. 

The most simple situation where (|1.2p is satisfied is when there exists a uniform 
gap between the eigenvalues: for all p G N*, E'p+i — Ep > Cq > 0. Note that in this 
case we have no = 0. This property is true in the following examples. 

- IfVciz) = a'^z^ + Vi{z), with ||Vi||l- < 2^/^?T^. Indeed, in this case the 
perturbation theory gives \Ep — {2p + 1)\/ + B'^\ < \\Vi\\l^. 

- If Vc{z) ~ al^l** as l^l — > +00, with s > 2. Indeed, in this case the Weyl 

2s 

asymptotics [19] gives Ep ~ Cp''+^ , so -Ep+i — Ep —> +00 as p — > +00. 

Let us now give a few indications on the Cauchy problem for (II. 1|) . (II. 2|) . (II. 3|) . 
This system benefits from two conservation laws, the mass and energy conservations: 

Vt>0, ||vi/-(t)||2, = ||.i,o||2^, £{^'it))=£{^o), (1.20) 

where the total energy of the wavefunction ^'^ is defined by 

Sm = ^Wdz^'Wh + ^WVVc^'Wh + ^\\{ed., + iBz)^'\\l, 

+ \\dy^lh + l\\Vv^^lh- (1.21) 

For fixed e > 0, the Cauchy theory for the Schrodinger-Poisson with a constant 
uniform magnetic field was solved in [TUdB] in the energy space. It is not difficult to 
adapt these proofs (see also the reference book [13]) to our case where an additional 
confinement potential is applied. The energy space in our situation is the set of 
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functions u such that £{u) is finite: 

= |u e L^{R^) : d,u G L^{R^), ^/VcU G L'^{M.^), dyU G L^{R^) 

and (d, + —]ueL\R^) 



This space seems to depend on e, which would not be convenient for our asymptotic 
analysis. In fact, it does not. Indeed, thanks to our assumption (|1.17p on the 
confinement potential, one has 



\ZU\\]^2 



< 



L2> 



so u G -B^ implies that zu G Lp' and thus dxU G L^. Hence one has 



VcU G 



and, on this space, we will use the following norm independent of e: 



Mil 



2 



+ ||(-A,.,,)V2^||i. + ||(//,)V2 



\u 



1 + 



L2 

^cU\\\2 +B'^\\zu\\\2, 



(1.22) 



where we used the selfadjointness and the positivity of —Ax,y and of the operator 
Hz defined by (jl.lQp . and where I denotes the identity operator. In this paper, 
we will assume that the initial datum ^'o in (|1.2p belongs to this space B^. Then, 
for all e > 0, the system p.ip . (I1.2|l . (|1.3p admits a unique global solution ^'^ G 
C'^([0, +00), B^). Our aim is to analyze the asymptotic behavior of ^'^ as e ^ 0. 

We are now in position to state our main results. Here and throughout this paper, 
we will use the notation 



Vn G l1(M), (n) = / u{z)dz. 



Let us introduce the limit system. First define the following coefficients 



\/p G N, 



Olr. 



1-E 



C^BzXpXq 



En 



(1.23) 



(1.24) 



where we recall that {Ep, Xp)pGN is the complete sequence of eigenvalues and eigen- 
functions of the operator Hz defined by p.l9l) . Then, we introduce the follow- 
ing infinite dimensional, nonlinear and coupled differential system on the functions 

(j)p{t,x,y): 

-ap dl^p - dl<^p + Wcl>p , 0p(t = 0) = (M/Q Xp) , (1-25) 



VpG N, 



id, 



Wp 



w 




(1.26) 



Note that the convolution in p.26|) holds on the variables (x, y) G M^. The equation 
p.26|l is nothing but the Poisson equation for a measure valued distribution of mass 
whose support is constrained to the plane z = 0: 



^■n\J + + z^ 



J2 



P5 



z=0 



z=Q 
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In order to compare with , we introduce the following functions: 

<^{t,x,y,z) = ^(t)p{t,x,y)xp{z), '^lpp{t,x,y, z) = ^^6''^'^^/''^ (j)p{t,x,y)xp{z). 

(1.27) 

Remark that "^app can be deduced from ^ through the application of the operator 
^itH^/e^ ^ unitary on B^: 

^ app ^ ^ • 

This explicit relation is the only dependency in e of the limit system (|1.25p . p.26|l . 
p.27|l . Our main result is the following theorem. 

Theorem 1.3. Assume that Vc satisfies Assumptions \l.l\ and \1.2\ and let ^'o £ -B^. 
For all £ G (0, 1], denote by £ C°([0, +oo),B^) the unique global solution of the 
initial system (jl.ip . (|1.2p . (|1.3p . Then the following holds true, 
(i) The limit system (ll.25p . (|1.26p . (|1.27ll admits a unique maximal solution ^%pp G 
C^{[Q,Tmax),B^), where Tmax G (0, +oo] is independent of e. IfTmax < +oo then 

ll*app(*> Ollfii ^ +00 as t^ Tmax- 

(a) For all T £ {0,Tfnax), we have 

lim 11^^ — II =0 

Comments on Theorem 11.31 

1. The cyclotron effective mass. Theorem [L3] thus states that, on all time intervals 
where the limit system (11.27(1 . (I1.25p . (I1.26P is well-posed, the solution ^''^ of the 
singularly perturbed system (II. 1|) . (|1.2j) . (jl.Sp is close to As expected, the 

dynamics in the y direction, ie parallel to the magnetic field, is not affected by the 
magnetic field, since the operator is still —dy. On the other hand, the situation 
is different in the direction x and the averaging of the cyclotron motion results 
in a multiplication of the operator —d'^ by the factor Op which only depends on 
Vc and B. The coefficient ^ plays in (|1.25p the role of an effective mass in the 
direction perpendicular to the magnetic field. We find that the effective mass in the 
Schrodinger equation for the mode p depends on the index p of this mode. We do 
not know whether these coefficients are positive for a general Vc- 

Notice that the effective mass could be predicted heuristically by the following 
argument. Denoting by kx, ky the wavevectors of the 2DEG in the plane the 
electron dispersion relation Ep{kx,ky) in the transversal subbands can be written 
from p.lj) by computing the eigenvalues of the operator 

^ ^ ^ +B'^z'^ + Vc{z) + 2eBzkx + e^kl + e^kl 



2 V dz2 



Since e is small, an approximation of Ep{kx, ky) can be computed thanks to pertur- 
bation theory, which gives the following parabolic band approximation: 

Epikx, ky) = ^-klY. i^BzXpX,)' +kl + kl + o(l). 

^ g^p ^I-EP 

We can read on this formula that the effective mass is 1 in the y direction and is 
according to (|1.24|) in the x direction. Note that the specific case of the harmonic 
potential is treated below (see comment 3). 

2. Conservation of the energy for the limit system. Let us write the conservation 
of the energy for the limit system. The total energy for this system can be splitted 
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into a confinement energy Sconfi^) and a transport energy Stri^) defined by 

£conf{<i>) = Y.EpUp\\l2, (1.28) 

peN 

^tri*^) = + ^\\9y(l>p\\h 

peN peN 

/ , /I L ===\M^,y)\^\M^',y')\^dxdydx'dy'. {1.29) 
^ ■^i^* ^tty'Ix - x'\^ + \y - y \ 

An interesting property is that these two quantities are separately conserved by the 
hmit system. If ^1^^ solves (fr25D . ([OS), (fT^Tll . then, for all t £ [0,r], we have 

£confi^lpp{t)) = Sconfi^lppm and Stri^lppit)) = ^^,(^'^^^(0)). (1.30) 

In particular, by summing up the two equalities in (|1.30|) . we obtain the following 
conservation property: 

£conf{Kpp{t)) + £tr{%j,p{t)) = Sconfi^ppiO)) + Stri^ppiO))- (1.31) 

Note that, in the general case, we do not know whether the energy defined by (ll.29p 
is the sum of nonnegative terms. This point is related to the fact that the well- 
posedness for t S [0, +00) of the Cauchy problem for the nonlinear system p.25|l . 
p.26|l is an open issue. Nevertheless, when the ctp are such that the energy is 
coercive on B^, ie when we have 

y^GB\ Co\mli<£conf{'^)+£tr{^)<Cimli+C2\m%i, (1.32) 

with a constant Cq > independent of e, then the maximal solution of p.25|l . (|1.26p 
is globally defined: Tmax = +00. 

Corollary 1.4 (Global in time convergence). Under the assumptions of Theorem 
assume moreover that there exists < a < a such that the coefficients ap 
defined by p.24|) satisfy the following condition: 

Vp G N, a<ap<a. (1.33) 

Then the system (ll.27p . (I1.25p . (I1.26P admits a unique global solution '^app ^ 
C°([0, +c)o), B^) and, for all T > 0, we have 

lim 11^^ — II = n 

imi||4' ^app||cO([o,T],Si) ^' 

where G C°([0, +00), S^) denotes the solution 0/(0), (fOll . 

The proof of this corollary is immediate and will not be detailed in this paper. 
Indeed, remarking that (|1.33|) implies (|1.32|) . we obtain that the solution ^^^^(t) of 
p.25|) . (|1.26|) satisfies the following uniform bound: 

ll^app(i)ll|i < C {£confmpp{t))+£trmppit))) = C (f^eon/ (^o) + ^tr (^o)) , 

where the quantity in the right-hand side is finite as soon as ^'o G S^. 

3. Case of harmonic confinement. In the special case of a harmonic confinement 
potential Vc{z) = a'^z'^, the eigenvalues and eigenfunctions of Hz = —d1 + {a?'+B'^)z'^ 
can be computed explicitely and one has 

Ep = {2p + l)Va^ + B^ xp{z) = (a^ + ^2)1/8 (^(^2 ^ ^2)1/4^^ ^ 
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where (Mp)peN are the normahzed Hermite functions defined e.g. in pi], B 8 and 
satisfying —u" + z'^Up = {2p + l)up. The properties of the Hermite functions give 



V2(p+1) ^ 



and one can compute expHcitely the coefficients 



,2(22XpXp+i)^ , D2(2^XpXp-i)^ _ 



Ep+i — Ep Ep — Ep^i a? + B"^ 

We thus recover here the coefficient found in subsection 11.31 in the simphfied situa- 
tion. Note that, in this case, condition (|1.33|) is satisfied and the convergence result 
holds on an arbitrary time interval. It is reasonable to conjecture that this condition 
(|1.33ll holds again when Vc{z) = a^z^ + Vi{z)^ where Vi is a small perturbation. 

^. Towards a more realistic model. Since we aim at describing the transport of 
electrons, which are fermions, our model should not be restricted to a pure quantum 
state. The following model describes the transport of an electron gas in a mixed 
quantum state and is more realistic: 

idt^ = ^ (-^2 + B^z^ + Vciz)) - -2iBzd^^ - J\^^y^ + , Vj, (1.34) 

^'^(0, X, y, z) = ^'j- o(x, y, z), Vj, (1.35) 

y^(^,x,.) = ^*p^ P^ = Y.^0f^ (1-36) 

3 

where Aj, the occupation factor of the state takes into account the statistics 
of the electron ensemble and is fixed once for all at the initial time. Note that the 
Schrodinger equations (ll.34p are only coupled through the selfconsistent Poisson 
potential. Therefore, we claim that our main Theorem 11.31 which has been given 
for the sake of simplicity in the case of pure quantum state, can be extended to 
this system p.34|) . p.35|) . p.36p . with appropriate assumptions on the initial data 
(*.,o) ■ 

Similarly, a given smooth external potential could be incorporated in the initial 
system. We also claim that our result can be easily adapted if we add in the right- 
hand side of (jl.ip a term of the form Fe2,.t(i, 2;, y, e-z)^^ (which is coherent with our 
scaling), and the result does not change qualitatively. 

1.5. Scheme of the proof. In this section, we sketch the main steps of the proof 
of the main theorem. 

Ste'p 1: a priori estimates. 

The first task is to obtain uniform in e a priori estimates for the solution of (ll.ljl . 
p.2p . p.3|l . which are of course crucial in the subsequent nonlinear analysis. Due to 
the presence of the singular and ^ terms in (II. ip . this task is not obvious here. In 
subsection 12.11 we introduce a well adapted functional framework: a Sobolev scale 
based on the operators —^x,y and Hz. More precisely, for all m S N, we introduce 
the Hilbert space 

B^ = [u: hill™ = h||22(^3) + ||(-A,.,j,)"^/\||2 2(^3) + ||//r/2^||i2(K3) < +00} . 

(1-37) 

In subsection 12.11 we give some equivalent norms which are easier to handle here. 
Then in subsection 12.21 we take advantage of this functional framework and derive 
some a priori estimates for (jl.ip , (11.21) , (|1.3p . 
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Step 2: the filtered system. 

In [HlEj, the asymptotics of NLS equations under the form 

1 



idtu" = H.u" - A^^yu" + J^(|n^^u^ (1.38) 



£ 



such as the Gross-Pitaevskii equation, was analyzed. In p.38|) . : M+ M is a 
given function and the nonhnearity depends locally on the density It appeared 
in [6] that a fruitful strategy is to filter out the oscillations in time induced by the 
term Hz, without projecting on the eigenmodes of Hz . Indeed, projecting (|1.38p 
on the Hilbert basis Xp leads to difficult problems of series summations and of small 
denominators in oscillating phases. Introducing the new unknown: 

i;^(t, X, z) = exp [itHz/e^) u'^it-, x, z), 

the filtered system associated to (|1.38p reads 



2 ^-itH,/s\s (J3g) 

where we used the fact that Hz, thus e**''^^, commutes with dx and dy. Then, the 
analysis of the limit e — > amounts to prove that it is possible to define an average 
of the nonlinearity in (|1.39|) with respect to the fast variable t/e^. 
Let us adapt this strategy to our problem. Introduce 

^%t,x,z) = exp [itHz/e^] ^'^(t,a;,z). 

One deduces from (jl.lj) . (|1.2p . p.3|) the following equation for 

where we introduced the nonlinear function 

{t,u) ^F{t,u) = e"^' * \e-'^^'uA e'^^^^u, (1.41) 



and where is still defined by (|1.4I) . 

Step 3: approximation by an intermediate system. 

Before performing the limit e — > in (|1.40l) . we remark that (|1.41l) can be approx- 
imated in order to get rid of the fast time variable t/e^ in the nonlinear term of 
(|1.40|) . By writing formally 

^ ^ +0(1), (1-42) 



we remark that 



L^.(|u|^)+o(l), 



where the symbole * denotes here a convolution in the {x,y) variables only, and 
where we used the fact that e^"^^^ is unitary on L^(R). Hence, inserting this Ansatz 
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in (fOTI) yields 

F{t,u) =e'^^ 




Denoting 




and introducing the solution of the following intermediate system: 

id^^e = (e**^^^/^'ze-**^^/^') - A,,j,r^ + Fq , (1.44) 
we expect that the solution ^''^ of (|1.40|) satisfies 

^^ = W^ + o{l). (1.45) 

Subsection 12.31 is devoted to the rigorous proof of this heuristics. We give sense to 
the o(l) in Lemma [2^ and we prove that the solutions of the two nonlinear equations 
(|1.40l) and (|1.44l) are close together and that (|1.45p holds true in the sense of the 
norm. This statement is given in Proposition 12.11 

Step ^: second order averaging of oscillating systems. 

Thanks to this Step 3, we can consider the simplest system (|1.44l) instead of (|1.40l) . 
We are now left with the analysis of the asymptotics of this intermediate system as 
e — > 0. Note that (11.44p is under the general form 

idtu=-f(^u{t)+g{u{t)) (1.46) 



with 

/(r) = -2Be'^"^ze-'^"Hd^ and g{u) = -A^^yU + Fo{u). 
At this point, a critical fact has to be noticed. Equations under the form 

idtu = f (J^^ u{t) + g{u{t)) (1.47) 

can be averaged when, due to some ergodicity property, one can give a sense to the 
time average 

/° = ^lim i / fir) dr. (1.48) 

Indeed, under rather general assumptions, the techniques of averaging of dynamical 
systems - see the reference book on the topic by Sanders and Verhulst [SZ]- enable 
to show that (|1.47|) is well approximated by the averaged equation 

idtu = f\{t)+g{u{t)). 

Yet, the oscillating term in (|1.46p . compared to the same term in (ll.47p . is multiplied 
by p Therefore, a necessary condition in order to perform the averaging of (I1.46P 
is that the average /° of / is zero. In our case, the integral kernel of the operator 
^itH^^^-itH^ ^ defined by 

Vu, e'^^^ze~'^"^u= [ G{t, z, z')u{z')dz' , 
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is given by 

In the last inequality, we used the fact that, by Assumption ll.il Vc is even. Indeed, 
this property implies that, for all p, {xpY is also even, thus {z{xpY) = 0. Conse- 
quently, since p q implies Ep ^ Eg, the kernel G{t,z,z') is a series of functions 
which all have a vanishing average in time. We thus expect that the operator- valued 
function /(r) has the same property: 

/°= hm i r fiT)dr = 0. 
T^+oo 1 Jq 

In such a situation, the theory of averaging has to be pushed to the second order 
[27] in order to obtain the limit of (|1.46p as e ^ 0. Section [3] is devoted to this 
question of second order averaging, which leads to the limit system p.25|) . p.26|l . 
The main result of this Section [3] is Proposition 13.21 

In the short last Section [H we prove our main Theorem [L3] by just gathering the 
results proved in the previous sections. 



2. The nonlinear analysis 

In this section, we obtain some a priori estimates uniform in e for the initial 
system (|1.2p . (11.31) and we prove that it can be approximated by an interme- 

diate system, where we regularize the initial data and where we replace the Poisson 
nonlinearity by its formal limit given in p. 431) . This intermediate system takes the 
form 

idt¥^ = irH^^^ - -2iBzd^¥^ - A^y¥^ + , (2.1) 

¥{0,x,y,z) =¥o{x,y,z), (2.2) 
W%t,x,z) = I * . (2.3) 

Notice that (|2.3p is nothing but the Poisson equation (|1.3p where we replace = 
\/ x'^ + y"^ + e^z"^ by = ^Jx'^ + y'^. Moreover, the initial datum ^'o in (12. 2p will 
be chosen as a regularization in of the initial datum ^q. Recall the definition 
(|1.37l) of the space B^. The main result of this section is the following proposition. 

Proposition 2.1 (Approximation of the initial system). Assume that Vc satis- 
fies Assumptions II. ij . 11.21 and that ^'o G B^ . For all e G (0,1], denote by 
'^'^ G C'^(M+,S^) the unique global solution of the initial system p.lj) . p.2p . (|1.3|) . 
Then the following holds true. 

(i) There exists a maximal positive time such that '^^ is bounded uniformly in e : 
the quantity 

ro:=supir>0: sup ||^'^||co([o,T],iJi) < +oo I . (2.4) 

I ££(0,1] I 



satisfies Tq G (0, -|-oo] . If Tq < +oo the 

lim sup||^'^||co([o,To],Bi) = +00. 
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(ii) For all T G (0, Tq), where Tq is defined by (|2.4|) . for all 6 > and for all 
integer m > 2, there exist ^'o € -B™ and es such that the following holds true. For 
all e G (0,55], the intermediate system (12. ip . (|2.2|) . (|2.3|) admits a unique solution 
^'^ G C'^([0, T], S™") satisfying the following uniform estimates: 

ye<es 11^^ -^^||co([o,T],Bi) < ^ (2.5) 

< Cdl^obOll^ob'". (2.6) 

Remark 2.2. It is a priori not excluded that Tq < +00. Indeed, although we are in 
a repulsive case, the energy conservation does not enable to obtain e-independant 
a priori estimates in (see the proof of Lemma \2M . This may be linked to 
the possible formation of caustics, as for the nonlinear Schrodinger equation in 
semiclassical regime, see e.g. 

2.1. Preliminaries. As we explained in subsection ll.51 our nonlinear analysis will 
deeply rely on the use of the functional spaces defined by p.37|) and adapted 
to the operators and —Ax^y. The following result was proved in [6j by using an 
appropriate Weyl-Hormander pseudodifferential calculus, inspired by [9l [22]: 

Lemma 2.3 ([6j). Under Assumption consider the Hilbert space B"^ defined 
by (I1.37P for m G N. Then the norm \\ ■ in (|1.37p is equivalent to the following 
norm: 

\W\\h"-(R«) + \\yc{z)"^^^u\\L2(M,3y (2.7) 

Moreover, for all u G B"^^^ , we have 

+ ||<9a;U||Bm + ||5yn||Bm + IIS^-uIIb™ + II \/K:'"II-B'" ^ ||^||b™+i- (2.8) 

The operator A^^y commutes with the rapidly oscillating operator e^**^^/"^^ and 
with the operator izdx- This will enable us to obtain uniform bounds for the solution 
of (jl.ip by simply applying Ax^y to this equation. Unfortunately, the operator 
does not satisfy this property. For this reason, we introduce the following operator: 

He = H^- 2ieBzdx - e^dl = -dl + Vc{z) + {iedx - Bzf . (2.9) 

This operator enables to define another norm equivalent to the B^ norm. The 
following lemma is proved in the Appendix lAl 

Lemma 2.4. The operator defined by (12.91) on L^(R^) with domain B^ is self- 
adjoint and nonnegative. There exists a constant Ci > such that, for all e G (0, 1] 
and for all u G B^ , we have 

^Ml^ < MhiR^) + \\{-Ax,y)'/MUM3) + \\H^^Mh(^R3) < Ci\\u\\l, . (2.10) 

Moreover, for all integer m > 2, there exists Em G (0, 1] such that, for all e G (0, Em], 
for all u G B"^ , we have 

\\\u\\l^ < ||n||i2(j,3) + ||(-A,,,)"^/2u||2,^^3^ + \\Hf\\\l,^^,^ < 2\\u\\l^ . (2.11) 

2.2. A priori estimates. In this subsection, we obtain a priori estimate uniform in 
e for the initial Schrodinger-Poisson model p.ip . (|1.3|1 and the intermediate model 
(|2.ip . (|2.2|) . (|2.3p . Remark first that these two models can be considered in a unified 
way. For all n G -B^ and for a G {0, 1}, denote 
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where the convolution holds on the three variables {x,y,z) G M^. Remark that for 
a = 0, this definition coincides with the definition (|1.43p . We shall consider for 
e G (0, 1] and a G {0, 1} the nonlinear equation 

idf^u' = ^Heu" - dlu" + FJv!") , (2.13) 

u%0, X, y, z) = uo{x, y, z), (2.14) 
where the operator Hgr was defined by (|2.9p . Note that for uq = ^'o and a = 1, 
(l2l3D . (I2lill is the initial system {TI]), ([Ol), dLS), and that for uq = and 
a = 0, (I2l3]l . (I2lill is the intermediate system ([21]), ([221), (l23|). Let us first 
state a technical lemma concerning the nonlinearities Fi and Fq, which is proved in 
Appendix [B1 

Lemma 2.5. There exists a constant C > such that, for all e G (0, 1], for a = 

or 1, we have 

yu,veB\ \\F^iu)-F^iv)\\s^<C{\\u\\% + \\v\\%)\\u-v\\Bi, (2.15) 

where F^ is defined by (|2.12|] . Moreover, for all m £ N* , there exists Cm > such 
that we have the tame estimate 

VeG(0,l],VaG {0,1}, Vugs'", ||F„(n)||^„ < C^||n|||i . (2.16) 

Now we are able to derive uniform a priori estimates for the solution of (|2.13l) . 
(1214)1 . 

Lemma 2.6. Let e G (0,1], a G {0,1} and uq £ . Then the solution of the 
equation ()2.13p . (|2.14|) exists and is unique in C^{[0, +oo), B^) and the following 
uniform in e estimates hold true. 

(i) For all M > 0, there exist T > 0, only depending on M and UnoUsi, such that, 
for all e G (0, 1], we have 

lk1lco([o,T],Bi) < (1 + M)\\uo\\b^ . (2.17) 

(a) Let m > 2 an integer and assume that uq G B"^. Then, for all T > 0, we have 
the estimate 

Ve G {0,em], ll«^llco([o,f],B™) - C\\uo\\b^ exp (cf ||u^||^o([o,r],Bi)) • (2-18) 
where £m > is as in Lemma \2M 

Proof. Step 1: the Cauchy problem and the conservation laws. For any given e > 0, 
the existence and uniqueness of a maximal solution G C^{[0,T), B^) can be 
obtained by standard techniques [13]. We leave this first part of the proof to the 
reader. This solution satisfies both L^ and energy conservation laws: 

Vt>0, ||n^(t)||i2 = 11^0111,2 and £a{u%t)) = £a{uo), (2.19) 

where the energy is defined by 

£a{u) = -^{HsU,u)l2 + ||9yU,||^2 + ^{Fa{u),u)L2 

= -2 11^2^11^2 + ^llv^^^lliz + ^\\{ed^ + iBz)u\\\2 + \\dyu\\\2 + ^(F„(ii), u)i2 . 

We recall that the operator Hfr is defined by (|2.9|) . These conservation laws show 
that the solution is global, ie that T = +oo. Unfortunately, due to the 
terms in this expression, one cannot use the energy conservation to get uniform in 
e estimates. Instead, we will directly write the equations satisfied by dxU'^,dyU^ or 
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{H^Y^'^u'^ and use the standard L^-estimates for these equations and the fact that 
the self-adjoint operators H^^ dx and dy commute together. 

Step 2: estimate. This yields 

idt{Vx,yu'){t) = ^He{Vx^yu') - dliVx^yu') + Vx^y {F^{u')) 

and 

idt {hII^u') it) = ^H,{Hl'^u') - dl{Hll^u') + Hl'^ {F^iu')) . 

Hence, 

\\u%t)\\L2 + \\Vx,yU%t)\\L2 + \\Hl/\%t)\\L2 < \\U0\\l2 + \\Vx,yUo\\L2 + \\HI/^Uo\\l2 

+C 1^ {\\Va:,yFa{u%smL2 + || i?,^/^^^ (u^ (s)) || ^2 ) ds 

and, for e S (0, 1], the equivalence of norms given in Lemma [2!4l yields 

\\u%t)\\Bi < C\\uo\\b^ +C f \\F^{u'{s))\\b^ ds 

Jo 

< C\\uo\\bi +C [ ds, (2.20) 

Jo 

where we used (I2.15P with u = to estimate Fa{u^{s)). Hence, by applying the 
Gronwall lemma to the integral inequality (I2.20p . we prove Item (i) of the Lemma. 

Step 3: estimate. Let T > 0, m > 2, uq € -B™ and let e G (0,e,„], where 
< em < 1 as in Lemma [2!4l Since the operators Hfr and Aj;^^ commute together, 

H^^'^u^ satifies the following equation: 

idt (h^^u') it) = ^H,{Hf^u') - dl{Hf^u') + Hf" {F^{u')) , 
thus, for all t e [0,T], 



L2 < WHf'u.WL^ + f WHf' (Faiu^s))) h2 ds, 
Jo 

<C\\uo\\Bn.+C [ \\FUu'{s))\\B^ds 
Jo 

< C'llnollfim + C||n'^||^o([o,T],Bi) y \\u'^ is)\\B^ ds, (2.21) 

where we used Lemma [2!4l and the tame estimate (I2.16p . Similarly, —A^^yU^ satisfies 
the following equation: 

idt{-Ax,yu'){t) = ^H,{-Ax,yu') - dl{-Ax,yu') - A^^y (^^(n^)) 

and, using the definition of B"^ p.37|l and (I2.16P yields: 

\\{-Ax,yrl^u^{t)U2 < ||(-A.,,)™/2^o||l^ + f ll(-A.,,)"/2 {F^iu^s))) ^2 ds, 

Jo 

< C'II'WoIIb'" + C'||'"'^|lc'0{[o,T],_Bi) j \\u^{s)\\b"^ ds. (2.22) 
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Therefore, by using again the equivalence of norms given by Lemma [231 and the L? 
conservation law in (|2.19p . we deduce from (|2.2ip and (|2.22|) that, for t < T, we 
have ^ 

Ik^Wlls'" < C'llnollsm + C||n^||^o([o^2.]^5i) J ||m''(s)||b'" ds, 
and the Gronwall lemma gives (|2.18|) . □ 

2.3. Proof of Proposition [2TTI In this subsection, we prove Proposition 12.11 ie 
we show that this solution can be uniformly approximated by a regular solution of 
the intermediate system. We first state a technical lemma on the Poisson kernels, 
which is proved in the Appendix O 

Lemma 2.7. There exists a constant C > such that, for all e G (0, 1], we have 

yueB\ ||Fi(u)-Fo(u)||^i < C7ei/3||n|||2, (2.23) 
where Fq and Fi are defined by (|2.12p . 
We are now ready to prove the main result of this section. 

Proof of Proposition 12.11 Let "ifo ^ B^, let an integer m > 2 be fixed, and define 
the regularized initial datum ^'o by 

¥o = {I- r?A,,,)-'"/' (/ + *o , (2.24) 

where r/ > is a small parameter that will be fixed further and where I denotes the 
identity operator. Denote by the solution of the initial system (jl.lj) . (|1.2|) . (|1.3p 
and by the solution of the intermediate system (|2.ip . (|2.2|) . (|2.3p with the initial 
datum (|2.24D . We shall estimate the difference 

Step 1: uniform hounds for . Let < e < 1. Prom Lemma [2?6l (i), we first deduce 
that there exists Ti > only depending on ||^'o||_Bi such that, for all e € (0, 1] 

ll^^llcO([0,Ti],Bi) < 2||^'o||bi- 

This implies that Tq defined by (|2.4|) satisfies Tq > Ti > 0. Clearly, if Tq < +oo, we 
have 

lim sup||^'^||co([o,To],Bi) = 

otherwise by reiterating the above procedure we could find a uniform bound on 
[0,r2] with T2 > To. 

Now we fix T G (0, Tq) and 5 > for the sequel of this proof. Definition (|2.4|) of 
Tq implies that 

\W\\co{[o,T],B^) < C(II*o||bi), independent of e G (0,1]. (2.25) 

Step 2: hounds for the initial datum ^q. First, we deduce from (|2.24p that 

(/ - A,.,, + H^fl^^Q = {I- r?A,,^)-'"/2(/ + riH,r^/\l - A,,, + H,)^/^^ , 
hence 

\\{I-K.,y + H,YI^¥4L^ 

< 11(1 - 77A,,j,)-™/2(/ + rjH,)-^'\l - A,,y + H,)'/HoU2 

< \\{I-A,,y + H,)'/^^o\\L^ 

where we used the fact that the operators (I — r/Aa;^^)""^/^ and (/ + r]Hz)~'^^'^ are 
bounded on L^, with bounds equal to 1. Therefore, using p.22|) . we obtain 

ll^bi < ll^obi, (2.26) 
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where we recall that the right-hand side is independent of e. 

Next, we get from (|2.24p the two following identities: for all integer £ < m, 

and 

Thus, from the bound 

we deduce that both operators {-Ax,yY^'^{I - 'r\^x,y)~^^'^ and Hi^'^{I + rfH^)"^/'^ 
are bounded on L^, with bounds equal to Crj"^/"^ , and thus 

<Cr''^\\^4B^, (2.27) 

where we recall the definition (I1.37P of the norms. 
Finally, we obtain also from (|2.24|) that 

(/-A,,j,+/7,)i/2(M/o-ro) = (/-(/- r?A,,,)-"/^(/ + r^H^r^'A {I-A,,^y+H-,f'H^ . 



Decompose v = {I — A^^y + Hz^^'^^o on the Hilbert basis (Xp)pGN of eigenmodes 
of Hz -. 

v{x, y, z) = ^ Vp{x, y) Xp{z) 



and denote by Vp{C), G ^ , the Fourier transform of Vp{x,y). By (|1.22l) . we have 

ii^o-^oii|i = 5] / (i-(i + r?ien-™/2(i + r?i?p)— 

Hence, using that 

Y] I \vp{0\^dC = W^oWm < +00 (2.28) 



and that 



e R', yp en, lim ( 1 - (1 + r/|^|2)-"'/2(l + r?K)-™/2 ) = 0, 

we deduce from Lebesgue's dominated convergence theorem and from the conver- 
gence of the series in (I2.28|) that 

lim ll^'o - ^ollfji = 0. (2.29) 

Step 3: uniform a priori estimates for . Consider 

Tr, ■■= sup{r e (0,r] : Ve G (0,1], ||^'^||co([o,T,],Bi) < 2||^1|co([o,T],Bi)}- (2.30) 

Note that, from (I2.26P and Lemma \2M (i), we know that G (0,7"] is well-defined. 
Then, from Lemma [23] (H), we deduce the following estimate: 

VeG(0,e,„], < m, ll*^llco{[o,T^],iJ«+i) < C' (|I^1IcO([o,t„],bi)) ll^olls^+i 

< C (ll^'^llcO([0,T],Bi)) 

<C7(||^'obi)||^o||Bm (2.31) 

where we used dOOl) and KOB . 
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Step 4-' estimate of the difference vl/^ — vj/"^. Using the notations defined in (|2.9|1 and 
(I2l2]l . and satisfy (I2l3]l . (l2lill with a = 1, uq = ^'o and a = 0, lio = ^ 
respectively. The Duhamel formulation of these equations read respectively 

i-t 

^0 

^ ft 

Jo 

Hence, for all t £ [0,T^] and e G (0,6^], 

ll^^(t) - ^^^(t)bi < 11^0 - ^ollfii + r \\Fi{^%s)) - Fi{¥{s))\\Bids 

Jo 

+ / \\Fi{¥{s))-Fo{¥{s))\\Bids 
Jo 

< 11^0 - ^ollsi +C [ W^'is) - ¥{s)\\Bids + 
Jo 

where we used (ITTbII . (l2:25]l . (lOnll . (l2:28]l and (lOTTl with £ = 1, coupled to (12:2711 . 
Here C denotes a generic constant depending only on T and ||^'o||_bi- Hence, by the 
Gronwall lemma, we get, for all t £ [0,T^], 

\\^'{t) - ¥^{t)\\Bi < [W^o - ^oIIbi + e^^. (2.32) 

Now, according to (|2.29l) . we fix r/ such that 

11*0 - ^ollsie < mm ( 2' 3II* IIcO([o,t],bi; 
and, in a second step, we fix G (0, £m] such that 

c4/V^/2e^^<min(^,i||^leO([o,T],Bi)) 
From (I2.32p . we deduce that 

Vt E [0,r^], Ve G {0,es], ||*^(t) - ^'{t)\\B^ < min (^6, ^||^''|bo([o,T],si) ) • 
Therefore, we have 



(2.33) 



< ^ll*'llco([o,T],Bi)- (2.34) 

We claim that = T. Indeed, if < T, then, applying again Lemma 12.61 at 
and using (12.341) enables to find r > such that, for all e G (0, 1), 

\\^^\\cO{[Tr,,T,,+T],B^) < 2||^'^||cO([o,T],Bi)> 

which, together with (|2.34p . contradicts the definition (|2.30p of T^. Finally, (|2.33p 
gives (|2.5I) and (I2.3ip with i = m — 1 gives (12. 6p . The proof of Proposition 12.11 is 
complete. □ 
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3. Second order averaging 

In this section, we focus on the intermediate system (12. ip . (12. 2|) . (12. 3|) as e goes 
to zero. As we explained in subsection 11.51 it is interesting to consider the filtered 
version of this equation. Let ^'o £ -B™ be a given initial data, let ^"^ be the 
corresponding solution of (|2.1I1 . (12.21) . (|2.3p and set 

•) = exp {itH,/e') ¥^{t, ■). (3.1) 
This function satisfies the system 

idt^e = {e^'^^l'\e-''"^l'^^ {id^^^) - A^,y$^ + Fq ($"(t)) , (3.2) 

$^(t = 0) = $0, 

where Fq is defined by p.43|l . The advantage of this intermediate system, compared 
to (|1.40p is that the nonlinearity Fo($^) has no dependence in the fast variable p-. 

We will analyze the filtered system (13. 2p in the framework of second order aver- 
aging of fast oscillating ODEs under the form (|1.46l) -see p7]-, that we adapt here 
to our context of nonlinear PDEs. Recall that (£'p)pgN, (Xp)pgN are the complete 
families of eigenvalues and eigenfunctions of the operator and denote by Hp the 
spectral projector on Xp'- 

v<i> G l2(m3), np$ = {^xp)xp- 

Introduce now the following unbounded operator on L^(M^): 

A, = -dlY,apIip with «P = 1-E ^T-F^' - (^-^^ 

p>0 q^p P 

With this notation, the limit system (|1.25p . p.26|) . (I1.27P can be rewritten in a more 
compact form as 

idt^ = Ao<^>-d^^ + Foi'i>), q;{t = 0) = ^o. (3.4) 

We state the main results of this section in the following two propositions. 

Proposition 3.1. Assume that Vc satisfies Assumptions \1A\ and 11.21 Then the 
following properties hold true. 

(i) The unbounded operator Aq defined by (|3.3I1 on L^(M^) with the domain 

D{Ao) ={$ G L\R^) : dl^apUp<^ G L'^{R^)} 

p>0 

is selfadjoint. Moreover, the operator Aq satisfies 

y£ > 0, Vn G 52no+4+f WAQuWsi < C\\u\\ B2no+4+e (3.5) 



where uq is as in Assumption \1.2[ 

(a) Let ^'o G . The limit system (|3.4|) admits a unique maximal solution ^ G 

C°{[0,Tmax),B^)- IfTmax < +00 then ||$(t)||Bi +0O as t Tmax- 

Proposition 3.2 (Averaging of the intermediate system). Assume that Vc satisfies 
Assumptions \lA \ and \1.2{ Then there exists an integer m >2 such that the following 
holds true. For $o e -B"", we consider the solution $^ G C°([0, +oo), fi"") ofJOl 
and the maximal solution $ G C^{[0,Tmax), B^) of the limit system with as 
initial data: 

_ ft _ 

$(t) = e-^*(^o-^«') *o - M e-^(*-^)(^o-^^)Fo(«>(s))ds. (3.6) 
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We assume that there exist T £ {0,Tmax), eo > such that 

M := sup ||'&^||co([o_T],B™) < +00. (3.7) 

ee(0,eo] 

Then we have 

l|5^-5|lcO([0,T],i3i) <£<^M, (3.8) 

where Cm is independent of e. 



3.1. Well-posedness of the limit system. In this section, we prove Proposition 
EH 



Step 1. Basic properties of the operator Aq. First, from Vc{z) > aP'z^ , we deduce 
that the pih. eigenvalue of Hz is larger than the pth. eigenvalue of the harmonic 
oscillator + (a^ + B'^)z^: 



Vp G N, Ep> ^/a? + B'^{2p + 1). (3.9) 
From Assumption 11.21 we deduce that the coefficients ap in (|3.3p satisfy 

< 1 + C{l+pr ^'^Bzxpx,? = 1 + C{l+pr\\BzXp\\l2 

q>0 

1 /o 

where we used (13. 9|) and that ||i?zXpllL2 < ■ Now, consider a nonnegative 
integer £ and u in 2?^"°+'^+^. Let no be defined as in Assumption [L2l and decompose 
u over the Xp family which is orthogonal in L'^. 

\\Aou\\lt =Yal\\dlnpu\\l, 

p>0 

< c j;i?2no+2||n^^||2^^^^ < c ||/7r+'npu||^,+, 

p>0 p>0 

^ ^ ^ ||np?/||^2nQ+4+« = ^^11^11^2710+4+* 

p>0 

where we used Lemma [2?3l This proves (|3.5p . 

Furthermore, by passing to the limit as ^ +00 in the identity 

N N 

D{A), ^ap(92np$,npvl/)i2 =YapiUp'^,dlUp^)L2, 
p=0 p=0 

we obtain that the operator ^0 is symmetric. Moreover, the equation Ao^ + i^ = f 
admits a solution $ G -^(^0) for all / E L^(M^). Indeed, the projection of this 
equation on Xp reads 

-apd^cpp + i4>p = fp 

and this elliptic equation can obviously be solved for all fp € L^(M^). Therefore, 
by the standard criterion for selfadjointness p6], the operator ^0 is selfadjoint. We 
have proved the first part of Proposition 13.11 



Step 2. Well-posedness and stability of the limit system. The operator ^0 being 
selfadjoint, the Stone theorem can be applied and the operator —iA^ generates a 
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unitary group of continuous operators e"*"^"* on and also on B^. The Duhamel 
formulation of (|3.4p reads 

ft 

= e-^*(^o-^y)^o - i / e^*(*^")(^"-^')Fo($(s))(is (3.10) 

(recall that Aq and 9^ commute together). Since, by (|2.15p . the application Fq 
is locally Lipschitz continuous on B^, it is easy to prove by a standard fixed point 
technique that (jS.lOjl admits a unique maximal solution <I> G (^"^([O, Tmax), B^)- The 
details are left to the reader. Note that, if Tmax < +oo, then ||<I>(t)||5i +oo as 
t — > Tmax- Item (ii) of Proposition 13.11 is proved. □ 

Remark 3.3. In fact, this strategy of proof by a fixed point mapping leads to a 
stability result. For all r/ > and for all T £ {0,Tmax), there exists 6ri^T > such 
that the following holds true. For all ^'o satisfying 

11*0 -*0bi < Sr,,T, 

the equation (13. 6|) 

_ _ ft _ 

$(t) = e-^*(^o-^«)^o - i / e-'(*-")(^«-^')Fo($(s))ds 

Jo 

admits a unique solution $ € C°([0, T],S^) and we have 

sup \mt)-m\\Bi <v- (3.11) 

te[o,T] 

3.2. Proof of Proposition 13.21 This subsection is devoted to the proof of Propo- 
sition 13.21 which relies on a reformulation of the Duhamel formula for (13.21) . 

Step 1: reformulation of the Duhamel formula. Introduce the following family of 
unbounded self-adjoint operators on L^(M'^) 

Vr G M, a(r) = -2Be'''"^ ze-^'^^Hd^ (3.12) 

with domain B^. Note that, from (I1.17P and Lemma [231 we deduce that, for all 

yueB"^, VtgM, ||a(T)'u||i2 < C||'u||b2. (3.13) 
The Duhamel representation of (13. 2|) reads 

ft ft 
¥Ht) = ¥o-^J^a (^) ¥^{s)ds - i (-A,,j,i^(s) + Fo ($^(s))) ds. (3.14) 

Introduce the primitive of a: 

\/ueB'^, VtGM, A{t)u= [ a{s)uds, (3.15) 

Jo 

which is well-defined as a Riemann integral, thanks to (I3.13|) . and is such that 

yuGB"^, VrGM, \\A{t)u\\l2 < Ct\\u\\b2. (3.16) 

Now, we notice that if G C°([0, T], 5^), then by ([321) we have that dt^^ £ 
C^{[0, T], i?^). Hence one can integrate by parts in the first integral of (|3.14l) and, if 
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m > 4, the following expression holds true for all t £ [0, T], in the sense of functions 
in C0([0,r],L2): 



^ ^ (^) « (^) ^(^)'^« - (^) ^(*) 
+e a(^) (-A,,j,$^(s) + Fo($^(s)))ds, 
where we used (13. 2|) to evaluate idt^^. Finally, the Duhamel formula (I3.14|) becomes 
¥{t) =¥o + J^A a ¥{s)ds - ieA 

+e / A (4) (-A,,j,$^-(s) + Fo($^(s)))ds 

- A,,j,$^(s) + Fo ) ds. (3.17) 



5iep 2: approximation of the Duhamel formula. Denote 

^{t) = W^{t) + ieA ¥^{t) 
and rewrite (|3.17p as follows: 



0^ / |A(^)a(^) +i9^)<D-(s)ds-i / [-d''^^'is) + Fo[^'{s)))ds. 











10 

In this step, we prove that 



sup \\^^{t)-^^it)\\Bi <eCM 

iG[0,T] 



and that 



(3.18) 



(3.19) 



$-(t) = ^o-i I ( ^o^^(s) - 5^«'^(s) + i^o ( ^^(s) ) + (s) ) ds, (3.20) 



with 



sup II/^IIbi < Cm- 

ie[o,T] 



(3.21) 



In order to prove this claim, we state two technical lemmas which are proved in the 
Appendix [P] so that the proof would be more readable. 

Lemma 3.4. Let Vc satisfy Assumptions II. ll and \1.2[ Then, for all integer £, the 
operator A{t) defined by (|3.15p satisfies 



te[o,T] 



A{-^] u{t) 



< C'||n||c-o([o,T],B2"o+^+8)! 

(3.22) 



where no is as in Assumption \1.2\ and C is independent of e. 
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Lemma 3.5. Let Vc satisfy Assumptions \l.l\ and \1.2[ Let T > and m = Auq + 17. 

Let u £ C°([0,r],S'") such that dtu G C70([0, T], S"^-^). Then we have, for all 
EG {0,1], 



sup 

t6[0,T] 



A 



+ id"^ u{s)ds + ^ y AQu{s)ds 



< C£^\\u\\ (3.23) 



where Aq, a and A are respectively defined by (|3.3p . (|3.12p and (|3.15p and where 
\\u\\ denotes shortly ||^i||cO([o,T],B™) + ll^t'f^ll(70([o,T],B™-2) • 

In order to apply these lemmas, we need some bounds for #^ and dt^^. Let us 
fix m = 4no + 17, where no is as in Assumption 11.21 and assume that we have the 
uniform estimate (|3.7|1 . By (|2.8p . we deduce that 



llco([o,T],B™-2) + 



C0([0,T],_B'"-2) 



Moreover, from (|2.16p . we deduce that 

l|-^o(^^)||cO([0,T],B™) < C'm- 

Hence, from ([321), (l3:24ll and ([3211), we get 

l|f^i^^llcO([0,T],B™-2) < 



< Cm. (3.24) 



(3.25) 



(3.26) 



Therefore, applying Lemmas [Ml and [33] and using JMI, ^2M, and 
yield 

yl f ^"j $^(t) < Cm, 



e sup 

tG[0,T] 



sup 

te[o,T] 



A -2 -A,,,<i>-(t)+Fo <i>^(t) 



(3.27) 
(3.28) 



and 



sup 

t£[Q,T] 



A 



<eCM, (3.29) 



where we used that m > 4no + 17, thus in particular m > 4no + 13 and m > 2no + ll. 
Hence, from (|3.27p . we deduce (|3.19|1 and 



Hi^' - ^')\\cOi[0,T],B^) + 



<eCM, (3.30) 

CO([0,T],Bi) 

where we also used the estimate (|2.15|1 . Moreover, from (|3.5|) and (|3.27p . we get 

\\Ao{¥ - $^)||co([o,T],Bi) < eCM. (3.31) 

Finally, inserting JI^H]), (lOl . (lOTTl in (l3T8]l yields (lOl with the esti- 

mate ([3?2T]) . 

5iep 5; a stability result for the limit system. First notice that (|3.20p implies that 
satisfies in the strong sense the equation 

idt^^ = - 5^$? + Fo($^) + e/^ $(t = 0) = 

which has the following mild formulation: 



(3.32) 
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Apply now Proposition 13.11 (ii) with ^'o as initial data: there exists a maximal 
solution # G C^{[Q,Tmax), B^) to the equation (|3.6I1 . Assume that T is such that 
< T < Tmax- Substracting jM]) to (HO^ leads, for all t < T, to 

< f Fo($?(5))-Fo($(s)) „/s + e||r|bo([o,T],Bi) 
./n i> 



ds 



where we used (|2.15p . (|3.21|) and ||'&^||cO([o,t],_b1) ^ ^m- Therefore, the Gronwall 
lemma gives the estimate (|3.8p and the proof of Proposition 13.21 is complete. □ 



4. Proof of the main theorem 

This section is devoted to the proof of the main Theorem 11.31 Remark that the 
statement (i) is. already proved in Proposition 13.11 Let us prove the statement (ii) 
of Theorem 11.31 

Let ^-0 G B^. Denote by W £ C^{[0, +oo), B^) the solution of JTI]), (fOll 
and let Tq G (0, +oo] be the maximal time given by Proposition 12.11 (i). We also 
introduce the maximal solution $ G C^{[0,Tmax), B^) of the limit system (|3.4|1 . 
given by Proposition 13.11 Pick T such that 

< T < mm{To,Tmax) 

and let 77 > 0. 

Since T < Tmax, according to Remark 13.31 one can define (^r^/s.T > such that 
the following holds true. For all ^'o satisfying 

11^0 - ^ollfii < \/3,T, 

the equation (13. 6p admits a unique solution ^' G C^{[0,T], B^) and we have (13. lip : 

sup ||^>(t)-$(t)bi <7?/3. 

i6[0,T] 

Next, we fix m > 2 according to Proposition 13.21 and 5 > by 

(5 = min Q,5^/3_T^ . (4.1) 

Since T < Tq, Proposition 12.11 (ii) enables to choose ^'o G B^ and es such that 
the corresponding solution of the intermediate system (|2.ip . (|2.2|) . (|2.3p satisfies 
(1231) and for all e < es: 

ll^'-^1bo([o,T],Bi) <'^<| (4.2) 

and is bounded in C^{[0,T], B"^) uniformly with respect to e. 
Now, we remark that by (|4.2p this initial data ^0 satisfies 

11^0 - ^ollfii < <^ < \/3,T- 

Hence, Remark 13.31 gives that the solution <1> of the equation (13. 6p satisfies 

~ 7] 
||<I> - ^||cO([0,T],Bi) < 

or, equivalently, 

||e I <P-e *llc0([0,T],Bi) < (4.3j 
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Moreover, the uniform bound of ^'^ in C^{[0,T], B"^) enables to apply Proposition 
13.21 which gives that the function ^f"^ satisfies 

ll^e _ e-^*^^/^'$|bo([o,T],Bi) < < (4.4) 

for e < es, where $ solves ([3?6l) . Finally, (1421) . ([43]) and (l44l) yield the existence 
of eo such that, for all e S (0, eo] we have 

ll^'-e-'*^^/"'<i>||co([o,T],Bi) (4.5) 

To conclude, it remains to remark that Tq > Tmax- Indeed, if Tq < T^ax, then we 
have, by Proposition 12.11 (i), 

lim sup||^'^||co([o,To],i3i) = +00, 
which implies by (14.5(1 that 

hin ||$(T)bi =+oo. 
This contradicts Tq < Tmax- The proof of Theorem 11.31 is complete. □ 



Appendix A. Proof of Lemma 12.41 

First, by integrating by parts and applying Cauchy-Schwarz, we obtain 

||-Bz92,.m||^2 = / B'^z'^ \dxu\'^dxdydz = / {B'^z'^u){—d'^u)dxdydz<\\u\ 



B2- 



Hence, the first properties stated in the Lemma are obvious from the definition 
(|2.9p . and we shall only detail the proof of the equivalence of norms. 



Step 1: the case m = 1. From the definition (|2.9p and the assumption (I1.17P on Vc , 
we deduce that 

WH^^Mh = i(-dz + Vc)u,u)l2 + \\{ed.,+iBz)u\\l, 
= {{—d'^ + Vc)u,u)]^2 +i?^||zu||^2 + e^||9i.ii||^2 — 2eBlm{zu,dxu)i2 
1 

> ^ii-dz+yc)u,u)L2 + i—+B^)\\zu\\l2+e^\\dM\l2-2£B\\zu\\L2\\dxu\\l2 

> ^((-^z +Vc)u,u)l2 + -j\\zu\\l2 + ^2 ^4^2 ^^ll^^^lli^ 

> C\\Hy\\\l2 + Ce^\\dxu\\l2. 

Conversely, from (|1.22p and (|2.9|) . we estimate directly 

{H,u,u)l2 < C'\\Hl/\\\l,+C'e^\\dxu\\l2. 

For all e G (0, 1], this yields the equivalence of norms (|2.10|) . 

For m > 2, we will proceed by induction. For the clarity of the proof, let us 
introduce two notations. For m S N, we denote by {Vm) the property 

(Vm)' there exists > such that, for all e £ (0, and for all u £ B"^, we have 

"'■|U,I|2 ^ IU,I|2 _i_ II A'"/2„||2 _L II H'™/2„||2 ^9lU,l|2 

-\\U\\Bm S |F|Il2(r3) + \\^X,y ^IIl2(R3) + \\J^£ U\\l'2(R3) S ^\\U\\sm, 



■',y 

and by (Qm) the property 
(Qm): there exists Cm > such that, for all u £ B"^ and e £ (0, 1], 

the operator Am = ^{H^ — H'^) satisfies \{AmU,u)]^2 \ < Cm||u||^ 
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Note that the lemma will proved if we show that (Vm) holds true for all m > 0. 
Note also that, up to a possible modification of the sequence (em)mGN, this sequence 
can be chosen nonincreasing. 

Step 2: (Qm) implies {Vm). Let m > be fixed. From (Qm), we deduce that 

\\HT''^u\\\2 = {H^u,u)l2 = {H'!^u,u)l2 +e{AmU,u)L2 

= \\H!p^^u\\l2 + e{A^u,u)L2 , 

thus 

\\HT^^u\\l2 - eCm\\u\\lm < \\H^/^u\\l2 < \\H!!'/^u\\l2 + eCm\\u\\ln. . (A.l) 
Setting 

1 



2 Cm, 

we deduce directly from (I1.37P and (jA.ip that, for e < ffm, 

''"IU,I|2 ^ IU,I|2 _i_IIaW2„||2 _l W U^/'^^,\\'^ ^9lU,l|2 

We have proved (Vm)- 

Step 3: proof of (Qm) for m = and 1. For m = 0, choose = and (Qo) is 
obvious. Let us prove (Qi). From (|2.9p . we have 

He = H,+ eAi, with Ai = -2iBzd^ - ed^. (A.2) 

For all u G -B^, we have 

\{Aiu,u)l2\ = I - 2iB{dxU,zu)i^2 + e\\dxu\\\2\ < C{\\zu\\l2 + \\dxu\\l2) < Ci||u|||i, 
where we applied Cauchy-Schwarz and Lemma [231 We have proved (Qi). 

Step 4-' proof of (Qm) for m > 2. We shall now proceed by induction. Let m > 2 
and assume that (Qm-2) and (Qm-i) hold true. Let us prove (Qm)- We compute 

= {H, + eAi)H^-^ {H, + eAi ) 

= H,H^-^H, + eAiH^-^ + eFf 

= + eH,Am^2H, + eAiH^-^ + sH^-^Ai 

where we have applied (Qm-2)- Hence, denoting 

A^ = H,Am-2H, + AiH^~^ + H^-^Ai , (A.3) 

we obtain 

H^ = H^ + eArn 
and, for all u € S™, we get from the definition (|A.3p that 

\{AmU,u)L2\ < \{H,Am-2H^u,u)L2 \ + 2\{H^-\, Aiu)\l2 , 

where we used that H^~^ and the operator Ai defined by (IA.2j) are selfadjoint. It 
remains to estimate the two terms in the right-hand side of this inequality. The 
first one can be estimated as follows: 

\{HzAm~2HzU,u)]^2\ = \{Am-2HzU,HzU)i2\ < Cm~2\\H zU\\\^.2 < Cm-2\\u\\%m , 
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where we used (Qm-2) and (I2.8p . The second one can be estimated as follows: 



K 



m—l 



u, {idx){—2Bzu + idxu)) 



L2 



He^ {idxu),He^ {-2Bzu + idxu) 



L2 



< 



m—l 



idxu) 



He ' {-2Bzu + idxu) 



L2 



< C||(?a;U|| £Jm-l llziill^Jm-l + C||9x^||fjm-1 
^ C 1 1 li] I 5 

where we used that commutes with dx, the Cauchy-Schwarz inequality, the 
property (Vm-i) and, at the last step, (12.81) . Therefore, we have proved that 

\{AmU,u)L2\ < Cm\\u\\%,n , 

which proves (Qm)- The proof of the lemma is complete. 



Appendix B. Proof of Lemma 12.51 

For readability, we introduce in this appendix the following notation: 



y{x, y, z) £ R\ Va G {0, 1}, Ve G (0, 1), r^(x, y, z) = + + ae^ z'^ . 

With this notation, for all u £ B^, and a E {0, 1}, the nonlinearity Fa{u) defined 
in (fZT2]l reads 



1 



In order to prove the estimates stated in Lemma 12.51 we prove the following 
technical lemma on the Poisson nonlinearity. 

Lemma B.l. The following estimates hold. 

(i) There exists a positive constant C that does not depend on e £ (0, 1] or a G {0, 1} 
such that 

yu,v £ H^{M.^), —*{uv) <C\\ u\\ t/l \\v 



(B.l) 



(a) There exists a positive constant C that does not depend one £ (0, 1] or a G {0, 1} 
such that, if D denotes a derivative with respect to x, y or z, 



D \ — * (uv' 



<C||n||^i||^;||;,i (B.2) 



r 3 r 00 



(Hi) For any integer k, let f3 = {Px,Py,/3z) G 6e a multiinteger of length = 
fix + fiy + (3z = k and let = dx'^d^^ dz" he the associated derivative. Then there 
exists a positive constant Ck depending only on k such that 



yu G h\ ( — 



* \u\ 



r 3 T 00 



Proof. Noting that, for all {x^y) G 



— * {uv) ) {x,y, ■) 



< 



M2 ^y{x - x'Y + {y - y'Y 



< Ck\\u\\Hi\\u\\H''- 



\uv{x',y', ■)\ 



(B.3) 

dx'dy', 
(B.4) 
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we only need estimates for the convolution with 



m 



Here, we refer the 



reader to Lemma B.l of |7j where it was shown that for any / E L'p{M?) n L^(E^) 
with 2 < p < oo, the following bound holds: 



^3.2 _|_ y2 



*/ 



<Cr, 



LP 



(B.5) 



where 9 = p/{2p — 2). Moreover, from Cauchy-Schwarz and Sobolev embeddings, 
we deduce that for all p G [1, +00), 



\uvix,y, Ollj 





< 







Hx,y,-)\\LHR) \\v{x,y,' 



IL2( 



LP 



Combined with (IB.4P and (jB.Sp . this proves Item (i). 

In order to prove Item (ii)^ consider a first order derivative D with respect to x, y 
or z and let n, ?; G //^(M'^). Usual properties of the convolution give 

D \ — * (™) ) = — * -C (uv) = — * (D(u)v + uDiv)) . 

Using (|B.4I1 combined with the generalized Young formula gives 



— * {D{u)v + uD{v)) 



' a 



< 



r 3 r 00 



^3,2 _|_ y2 



* \\D{u)v + uD{v)\\^^ 



< C 



\D{u)v + uD{v) 



I LI 



r6/5 



(B.6) 



since the function x 



belongs to L 



2 (Ta2\ 



noting that, thanks to Sobolev embeddings. 



We end the proof of Item (ii) 



D{u)v + uD{v)\\^, < C\\D{u)\\l2\\v\\li^^lI+C\\D{v)\\l2\\u\\li^^lI 

< C||n||/^i ||f . 

In order to prove Item (in), we follow the same lines with derivatives of higher 
orders. Consider the derivative -D^ where /3 = {Px,f3y,Pz) G is a multiinteger of 
length \(3\ = Px + Py + (3z = k. Usual properties of the convolution gives 



Dn-*{\u\') 



L3 



L3 L°° 



Again, using (IB.4P combined to the generalized Young's formula lead to 

1 





< 


' a 


T 3 r 00 

^x,y-^z 



LU 



< c 



r6/5rl 

^x.y 



{B.l) 



We now write 

dI^{uu) = Cf3>D'^\u)D'^~'^'(u), 

I3'<I3 

where the sum is over the set of multiintegers /?' = (/3^, such that < 
Py < Py and P'^ < Pz- Thus, combining (|B.7I) with Sobolev embeddings gives as 
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above 



r 3 T oo 



< \\D^'i^)\\L4u\\Ll^Ll 

\/3'\=k 

\l3'\=i,0<£<k 

< C||n||j:^fc ||n||j|/i + C ^ ||'u||_H-*+i II^^IIh*-^- 



We conclude noting that, by interpolation, for all ^ < A; — 1, 



(B.8) 

□ 



Proof of LemmalKR We first prove (|2.15|1 . In that view, let us fix u and v in B^, 
£ G (0, 1) and a G {0, 1}, and note that 

1 



F^{u)-F^{v) 



47rr^ 



* + \v\){\u\ 



;|)])^+(J^*(H^))(n-.) (B.9) 



Hence (|2.15p is a straightforward consequence of the following claim. There exists 
a positive constant C such that for all ui,U2 and M3 G 

^ ' ' <C\\ui\\bi\\u2\\b4u3\\b^. (B.IO) 



* {UIU2) U3 



Proof of the claim (IB.lOp . According to Lemma 12.31 we have 



— * {U1U2) U3 



< C 



\ * {U1U2) ) U3 



-c 



( ^ * {U1U2) ) U3 



L2 

(B.ll) 



First, applying (IB.lj) and then Lemma [2^ 



— * (^1^2) lis 



< 



L2 



— * {UIU2) 



L2 < C\\ui\\bi\\u2\\b^ 



Similarly, we have 



< 



L2 



4r * {U1U2) 



< C||iii||Bi||n2| 



(B.12) 



(B.13) 



Moreover, if D denotes any differential operator of order 1, 



D 



1 



1 



* {U1U2) \u3 \ = D [—* (U1U2) ns + — * (U1U2) D{u3). (B.14) 



Applying the Holder inequality leads to 
1 



D 



* {U1U2) ns 



< 



L2 



* (^1^2) 



< C||ni||5i||n2||Bi 
where we used (IB.2j) . Finally, using (IB.lj) . 



(B.15) 



* {U1U2) D{U3) 



< 



L2 



* (U1U2) 



11^(^3)11^2 



<C||^xibi||n2bi||n3|| 



(B.16) 
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We d educe the claim ([RTn]l by combining (IbUI) with (IbH) . (IrTsI) . (IbUI) . (IrTKIi 
and (iRlHll . 

In order to prove (|2.16|1 . consider a positive integer m and fix n G i?™. According 
to Lemma [231 we only need to estimate and ||V"c™''^-Fa(ii)||i2. In that 

view, we readily have 

1 



< 



L2 



1/1 ' a 



< C'll^ll|i(R3)lhl|B' 



(B.17) 



where we applied (IB.lj) and Lemma [231 

Now, let denote any derivative of length m and write 



Hence, 



u 



L2 



l<|/3'|<m 



/3'</3 

1/3' I =m 



1 



1 



* \u\ 



\\D^-^'iu)\U,M (B.18) 



r 3 r oo 



< C\\ lill^l ll^ll/f™ + C ^ ^ Ijlill Jiff ll^lli^m-^+l 



where we applied (jB.ip . (IB.3|) and Sobolev embeddings. Using the interpolation 
estimate (|B.8|1 gives 

p™F„('u)||^2(K3) < C\\u\\]j4u\\Hm < C||^/|||i||n||B™ . (B.19) 

We conclude the proof of (I2.16P combining (|B.17p and (IB.19|) . This ends the proof 
of Lemma 12.51 □ 



Appendix C. Proof of Lemma 12.71 

In this section, we set for simplicity X = {x,y) G M?. In order to prove estimate 
(|2.23|) . we first study the difference between both convolution kernels. 
First Step: Difference between the convolution kernels 
Let u,v be two functions of B^. Denote 

^' = L L [ ^\x-xv^.H.-.r - \^\) 

(C.l) 

We split the integral as follows: 

5{u,v){X,z)=5+{u,v){X,z)+5-{u,v){X,z)= I / + 

Jx'en+ Jz'm JX'en- Jz'm 

where 

n+ = {X' £ \x - x'\ > e}, n- = {X' G \X - X'\ < e}. 

For all r/, /i G M, and X' ^ X, we have 

1 1 r'' 



y/\X - X>\^ + e^T]^ ^\X -X'\^ + e^fi^ Je^ (\X-X'\^ + C 



/ ^ Tj^d^ (C.2) 

Jeu f|X-X'|2+f2A3/2 
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and 



< 



Besides, a simple study gives 



< 



+ - 3^/3 



(C.3) 



(C.4) 



Equation (IC.2j) . combined with (IC.3j) and (|C.4p allows us to claim that for all 

0e (0,1), 



1 



1 



Now, applying (|C.5p with rj = z — z' , fi = z' and ^ = 3/8 leads to 
^1 1 ' 



\X -X' 




V|X-X'|2+e2(^_^/)2 y/|X-X'|2+e2, 



/2 



< Ce3/8|^|3/8 



1 



|X-X'|1V8 



|n(x',.Hx',-)llarf^' 



where we used the Holder inequality and Sobolev embeddings. Similarly, applying 
([OS]) with r/ = z' , ;u = and 6* = 3/4 leads to 



\u{X',z')\\v{X',z')\dz'dX' 



In+JR \ ^\X -X'\-^ + e^z'^ \X-X' 



z'f/\iX',-)viX' 



o+ \X-X>\V^ 
< Ce^/^\\u\\B2\\v 



dX' 



T2 



We have proved that 

v){X, z)\ < Ce^l\l + W'')\\u\\b. \\v\\b^ 
Consider now 5~ . Using (|C.2p again leads to 



(C.6) 



\5-{u,v){X,z)\ < 




(|X-X'|2+e2)3/2 



u{X',z')\\v{X\z')\didz'dX' 



Moreover, a simple computation gives 



(C.7) 



di 



(|X-X'|2+ ^2)3/2 \X-X>\' 



Hence, (|C.7p gives 

\5-[u,v){X,z)\ < C 




1 



-\u{X',z')\\v{X',z')\dz'dX' 



\X-X' 



< Ce'/^u\\Bi\\v\\Bi. 



(C.8) 
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Combining (|C.6j) and (IC.8j) allows to conclude that 



\6{u,v){X,z)\ <Ce'/^{l + ./V^)\\u\\B2\\v\\Bi, (C.9) 
where we have used z^^^ < C(l + y^VAz)), deduced from (fTTTl) . 

Step 2: Difference between the nonlinearities. In order to prove Lemma [2. 71 we need 
to estimate the following quantity in B^: 

Fi{u) - Fo{u) = 5{u,u) u, (C.IO) 

where u G is given. According to Lemma \2M we have 

\\Fi{u) - Fo{u)\\bi < C\\^c{Fi{u) - Fo{u))\\l2 + C\\Fi{u) - Fo{u)\\h^. 

First, we deduce from (IC.9|) that 

\\{l + y^c)6iu,u)u\\L2 <Ce^/^\\il + Vc)u\\L4u\\B4v\\B^ < Ce^/^Ml^, (C.ll) 

where we used Lemma 12.31 Let now D denote a first order derivative with respect 
to X, y 01 z. We clearly have 

\\D{Fi{u)-Fo{u))\\l2 



< 



— — I * {D{u)u + uD{u)) u 



+ \\5{u,u)D{u)\\l2 

L2 



< 2\5(u,D{u))u\\l2 + \\5iu,u)D{u)\\L2. (C.12) 
According to (|C.9|) . we have 

\\6{u, D{u))u\\l2 < Cei/3||(i + v^)u||i2 ||Mb2 < Ce^/^Ml^ (C.13) 

and 

< Ce^/mi + v^)Z)(n)||i2||nb2||nbi < Ce^/''\\u\\l„ (C.14) 

where we used again Lemma EHl Combining SCHM . (ICTTTIl . (IClTll . (IClTSll and 
(ICTTill gives ([223]). The proof of Lemma O is complete. □ 

Appendix D. Proof of the technical Lemmas 13.41 and 13.51 

Let us develop the operators a and A defined by (|3.12p and (I3.15P on the eigenbasis 
Xp- We have 

a{T)u = -Y.Yl e'^^^^'^'^apg id^Ug Xp 

p>0 q>0 

where we have introduced the coefficients 

ttpq = {2BzXpXq) ■ (D-1) 
Recall that, by Assumption 11.11 the potential Vc is even, so for all p, the function 
(Xp(^))^ is even. Therefore, we have 

Vp G N, app = {2Bzxl) = 0, 

thus 

a{r)u = -Y.Yl ^'"^''"''''''^(^P, ^dxUq Xp ■ (D.2) 

p>0 qj^p 

Let us now integrate this formula in order to compute the operator A defined by 
(I3J5I1 : 

A{T)u = i'^^ — ^ _^ apqid^UqXp- (D-3) 

p>0 qi^p 
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Before proving Lemmas 13.41 and f3.5l let us give a useful estimate on coefficients apq. 
For all p G N, g e N, A; G N we have 



^(fe+l)/2 
\0'pq\ ^ C' 



k/2 



(D.4) 



Indeed, we have 



2B 



rjk/2 

Hz' Xp,ZXq 



L2 



L2 



2B \^Xp,^-z 

<2B\\H^/\zx,)\\l^ 
< 2S|kx<?llBfc 



where we applied Lemma [2.31 



Proof of Lemma I3.4L Let no be as in Assumption 11.21 let £ G N and u G 

C0([0,r],52"o+8+^). Denoting 



— i'^Xp) ) /^p — ll'"pXpllc'0([0,T],S2no+8+^)' 



we have 



u 



CO([0,T],B2no+8+«) 



^fJ-l< +00. 



(D.5) 
(D.6) 



p>0 



From (Id31) . we obtain 



p>0 ijT^p 



no 



^pg 



ll'"gXpllcO([0,T],B«+i) 



where we used Assumption 11.21 Besides, applying Lemma [231 gives 

1 

ll'"gllcO([0,T],/i'»(R2)) 



E, 



< C- 



no+4+{£-s)/2 



,/2 



no+i+e/2 

Eg 

for all s < 2no + 8 + £. Hence, from the definition (|1.37p . we get 

lkgXpllcO([0,T],B<+i) < C!E^^'^'^^^'^\\Ug\\cO(^[o,T],L^{R'2)) + C'll'"gllcO([0,T]H<+i(IR2)) 

p{€+l)/2 p(^+l)/2 
^ ^P 

and, by using (ID^ and dSj]), 

(1+9)"" kpgl lkgXpllcO{[0,r],B'^+i) < 
Therefore, 



(D.7) 



I I 
p2 \^Pl\ 



E. 



(e+5)/2 ^ ^{e+i)/2 

TP2 



j-,no+4+£/2 

Eg 



- P2 TP 



El Eg 



1^ I 72 ) "(*)llcO([0,T],BO < 



\ 3/2 / 1/2 

by Cauchy-Schwarz. To conclude, it suffices to use (13.91) and (ID.6P : the series 
converge and we obtain the desired estimate (|3.22l) . 
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Proof of Lemma \EM Let m = Auq + 17 and let u G C°([0,r],S™) such that 
dtu e C°([0,r],S"^-2). Denoting now 

Up = {uXp) , l^p = \WpXp\\cO{[0,T],B^) + ll^t^pXpllc'0([0,T],B™-2)5 (D-8) 

we have 

ll'"llc'0([0,T],B'") + \\^M\cO([0,T],B^-2) ='^^p < +°°- (^-9) 

p>0 

Applying Lemma 12.31 as above yields 

4™"'^^'|I«PIIC0([0,T],//=(M2)) +i^i^-'~'^/'ll^^tnp|bo([o,T],H=(R2)) < CUp (D.IO) 

for all s < m. By composing the expressions (ID.3j) and (ID.2j) for and a, we obtain 

QiT{Ep-Eq) _ 

A{T)a{T)u =iY.Y.Y. ^^-r^^ e'^'^'^-^-^apqaqn dlun Xp 



Ep — Eg 
I _ gi^iEq-Ep) 



P^O QJ^P "7^9 ^ 



= ^ — ^ _^ — M"^ ^l^p 

p>0 q^p P 9 

P>0 gT^p n ^ g P 9 

Now, remark that, by (|1.24|) and (|D.lj) . we have for all p G N the identity 



g^p " 



Therefore we get, using the definition (|3.3 
(^(r)a(r) + id'^) u = — Mqu 



P>0 QT^p ^ 



(flp?)' 



P>0 5^P n^g P 9 

nytp 



and, integrating, 

y ^-^^ + a n(s)(is + i J AQu{s)ds 

= - EE#^XP fe-^^^-^^y^''dlup{s)ds (D.ll) 

In order to estimate the right-hand side of this identity, we claim that, for all p G N, 
p G N and A 7^ 0, we have 



Xp{z) l\'^'/'"dXis,x,y) ds 



CO([0,T],Bi) 1^1 

where Ct only depends on T and Un is defined by (jD.Sp . This claim is proved below. 
As a consequence, we can estimate (ID. lip as follows: 
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''^ r / s\ / S 

A 



£2 ; - V .2 



+ id^^ u{s)ds + W AQu{s)ds 



CO([0,T],Bi) 



^ 2^2^ |P _ P |2 ^{m-5)/2^P 

p>0 qjLp IP 91 ii/p 



p>0 (j^p „ ^ g 
ra T^p 



pl/2 pl/2 



2^2^ R2 (m-5)/2 
p>0 q>0 9 -'^P 



E,no+ll/2 



I V I ; p2p"0+2 (m-4)/2 " 

p>0(j>On>0 -^p-^? -C/n 



+ n 



< (^^2 1 _L rp2 ^ 1 

-^^ Z^Z^ (1 + ^2) 1+ 3^^^ Z^Z^Z^(l+p2Wi+ 2)1 

p>0 g>0 ^ ^ ' ^ p>0 g>0 n>0 ^ ' f > \ ' ) 

where we used Assumption II. 21 (ID.4p . (|3.9p and recall that m = 4no + 17. Hence we 
deduce (|3.23p by using Cauchy-Schwarz and (lD.9j) . It remains to prove the claim. 

Proof of the claim (lDl2l) . Let 

i-t 

v{t,x,y,z) =xp{z) e'^'^'' dlug{s,x,y) ds, (D.13) 
Jo 

for p G N, g G N and A 7^ 0. An integration by parts in (|D.13h yields 

v{t, X, y, z) = ijXp e'^'/'' d^dtUgis, x, y) ds + e^^*/"' dlug{t, x, y) - dlug{0, x, y)^ . 
Hence, by using (|D.10j) . we obtain 



bllcO([0,T],Bi) < Ct 



+ eT 



|A| 

where Ct only depends on T. This concludes the proof of (|D.12j) . 

The proof of Lemma [3.51 is complete. □ 
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